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Electronic Mechanism for the Coexistence of Ferroelectricity and Ferromagnetism 
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We study the strong coupling limit of a two-band Hubbard Hamiltonian that also includes an inter-orbital 
on-site repulsive interaction Uab- When the two bands have opposite parity and are quarter filled, we prove that 
the ground state is simultaneously fenomagnetic and ferroelectric for infinite intra-orbital Coulomb interactions 
Uaa and Ubb- We also show that this coexistence leads to a singular magnetoelectric effect. 

PACS numbers: 71.27.+a, 71.28.+d, 77.80.-e 



Mtwduction. The interplay between order parameters of 
different nature opens the door for designing new multifunc- 
tional devices whose properties can be manipulated with more 
than one physical field. For instance, the spin and orbital elec- 
tronic degrees of freedom can order individually or simultane- 
ously producing different phases. In particular, orbital order- 
ing can produce symmetry-breaking states like orbital mag- 
netism, ferroelectricity (FE), quadrupolar electric or magnetic 
ordering, and other multipolar orderings. The magnetoelec- 
tric multiferroics, such as i?(Fe,Mn)03 fl H H 0, H and 
i?Mn205 0.01, are real examples of materials that combine 
distinct useful properties within a single system. Recent stud- 
ies of these multiferroic materials have revived interest in the 
magnetoelectric effect, i.e., the induction of polarization with 
an applied magnetic field and magnetization with an applied 
electric field. To date, materials that exhibit both ferromag- 
netism (FM) and FE are rare 0] because the transition metal 
ion of typical perovskite ferroelectric s is in a nonmagnetic 
electronic configuration. Therefore, it is essential to ex- 
plore alternative routes to the coexistence of the FM and FE 

mM. 

Different mechanisms for FE involving electronic degrees 
of freedom have been proposed. There are those in which FE 
results from bond ordered states induced either by electron- 
phonon coupling (Peierls insta bility ) or by pure electron- 
electron Coulomb interactions 11211 . In these cases, the ferro- 
electric state is clearly nonmagnetic due to the singlet nature 
of the covalent bonds. In contrast, considering a system of 
interacting spinless fermions with two atomic orbitals of op- 
posite inversion symmetry (say the d and / orbitals), Porten- 
gen et al. Ql predicted that permanent electric dipoles are 
induced by spontaneous d — f hybridization when particle- 
hole pairs (excitons) undergo a Bose-Einstein condensation. 
This result was confirmed in the strong coupling limit of an 
extended Falicov-Kimball spinless fermion model where both 
bands are dispersive It was also confirmed numerically 
in the intermediate couplin g re gime by using a constrained 
path Monte Carlo approach 11511 . 

The critical question that emerges is: How do FE and 
magnetism interplay when real electrons, instead of spinless 
fermions, are considered? In this Letter, we answer this ques- 
tion, proving that the mechanism proposed by Portengen et 
al. I13I1 can coexist with magnetically ordered states. In this 



case, the single electron occupying the effective (say d — f 
hybridized) orbital simultaneously provides an electric and a 
magnetic dipole moment, and the Coulomb repulsion is suffi- 
cient to generate a strong coupling between both of them. 

We start from a two-band Hubbard Hamiltonian that in- 
cludes an inter-orbital on-site repulsive interaction Uab- Like 
in the spinless fermion case, this interaction provides the 
"glue" for the formation of excitons. At quarter filling and 
in the strong coupling limit, we map the low energy spectrum 
of the two-band Hubbard model, H, into an effective spin- 
pseudospin Hamiltonian, HcS, where the pseudospin repre- 
sents the orbital degree of freedom. We prove that in the limit 
of large intra-orbital repulsive interactions Uaa: Ubb — 00, 
Hcs has a ferromagnetic ground state that can be partially or 
fully saturated. By combinin g th is result with the previous 
analysis for spinless fermions Il4ll5[ll6ll . we show that FM 
and FE coexist, and a divergent magnetoelectric is demon- 
strated using the SO(4) symmetry of H. Our conclusions are 
reinforced by a semi-classical and a numerical computation of 
the zero temperature (T = 0) phase diagram of H that goes 
beyond the limiting case Uaai Ubb 00. 

Hamiltonian. We consider a two-band Hubbard model 
with a local inter-band Coulomb interaction Uab on a D- 
dimensional hypercubic lattice 

^ ~ ^i^t'' (/iLcr/i+e,,;>'cr + /i+e^i^'o-Zii^fj) 

i.,ri.,v,iy' ,a 

+ ^ Uabn1n\ + ^ C/^^n|'|n['| 



where r] = {x,y,z,...}, v = {a, 6}, = fl^fi^^, 
n'( = X!(T ^kr and rii — Since the two orbitals, a 

and b, have opposite parity under spatial inversion, the inter- 
band hybridization term must be odd under this operation: 
tab ~ —tba- In addition, the intra-band hoppings taa and 
tbb will have opposite signs in general. The local spin and 
pseudospin operators are given by the expressions: 

act' viy' 

(2) 

where ct^ are the Pauli matrices with /i = {x,y,z}- The 
total spin and pseudospin per site are: = X^jy ^iu ^'^'^ 
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rt^ = r/^. The pseudospin component is proportional 
to the on-site hybridization. Since the the orbitals a and b have 
opposite parity, the local electric dipole moment is pi = AtTj^, 
where /jt is the dipole matrix element between the a and h or- 
bitals Clfil. 



Symmetry is a useful concept for describing the coexis- 
tence of different order parameters. 101 Foriaa = ±^66, tab = 

ha = 0, Ea - = Cfc - Ubb/% and Uaa = Uab = Ubb, H 

is invariant under a U(1)(8)SU(4) symmetry group. The U(l) 
symmetry corresponds to the conservation of the total num- 
ber of particles. The generators of the SU(4) symmetry group 
are the three components of the total spin = J^i and 
pseudospin = J2i'''t P^^^ ^^^^ operators: 



E 



f- ,f; 



(3) 



The total spin is conserved for any set of parameters. If we just 
impose the condition tab = tba = 0, the symmetry group of 
H is reduced to the the subgroup U( 1 ) x U( 1 ) x SO(4). The six 
generators of the SO(4) group are the three components of to- 
tal spin st^ and the three operators tt^^. This symmetry arises 
from separate total spin and charge conservation of each band, 
as the two bands are only coupled by the Coulomb interaction 
Uab- The symmetry operators ± tt^^ are the generators of 
global spin rotations on each individual band (with the + sign 
for the a band and the — sisn for the b band). 



Strong Coupling Limit. We will consider from now on the 
quarter filled case Ua + rib ~ 1, where Ua and Ub are the par- 
ticle densities of the bands a and b. When Uaa, Ubb, Uab ^ 
\t,^^'\, the low energy spectrum of H can be mapped to an ef- 
fective spin-pseudospin Hamiltonian by means of a canonical 
transformation that eliminates the linear terms in t,.,,i : 



where H-^^^ 



J XZ 



Jn = 



Ubb 

8{tbbtaa 



i and 
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(4) 



tbb 1 

Thb^' 

In this limit, because 



Uab 

and = ea - £(, + {Ubb - Uaa)/"^ 
the double-occupancy is forbidden in the low energy Hilbert 
space of Hcs, both Si and Ti belong to the S* = 1/2 repre- 
sentation of the su(2) algebra. The first two terms of iJoff 
couple the spin and the orbital degrees of freedom. As usual, 
the Heisenberg antiferromagnetic interaction, H^^^, is a direct 
consequence of the Pauli exclusion principle. On the other 
hand, the anisotropic Heisenberg-like pseudospin-pseudospin 
interaction reflects the competition between an excitonic crys- 
tallization or staggered orbital ordering (SOO) induced by the 
Ising term, and a Bose-Einstein condensation of excitons in- 
duced by the XF-term The first term of i^cff shows 
explicitly that the amplitude of the excitonic kinetic energy (or 
Xy-pseudospin) term gets maximized when the excitons are 
in a fully polarized ferromagnetic spin state. However, anti- 
ferromagnetism (AF) is clearly favored by the second term. 

Large Uaa, Ubb limit. We will first prove that there are par- 
tially and fully polarized ferromagnetic ground states of iJcff 
in the limit of Uaa, Ubb oo and tab = 0, and that the total 
spin or magnetization st take the values < st < N/2. 
After proving this result, we will show that these ferromag- 
netic solutions are also ferroelectric for i?^^ < Bz < B'^ 
and, using the SO(4) symmetry, we will derive an exact ex- 
pression for ground state electric polarization as a function of 
the magnetization. 

Since Jo and Ji vanish in this limit, H^s is reduced to: 



cff 



E 

(iJ> 



[J^f r; + J^(t?t? + T^T^)H^^ + Bz ^ rf 



where the angular brackets indicate that the sum is over 
nearest-neighbors, J± ^ ~ Jy, and iJ/j = Si • sj + i. 
To prove our statement we will use a basis of eigenstates of 
the local operators t-^ and sf: {|rf ...t^) (g) |sf ...s^)} where 
N is the total number of sites. The off-diagonal matrix ele- 
ments of HcS 



(r-...Tf|ffeff|Tr...r 



N/ 



2J^5](r-...Tf|i^|^J|r^..r^)(.-....f|i^i^|.^..s^) (5) 

(iJ) 
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are non-positive because J± < 0, and the matrix elements of 
Hfy — T?T^ + and Hf^^ are explicitly non-negative. Ac- 
cording to the generalized Perron's theorem (see for instance 
Ref. LJ), there is one ground state of iJcff , 



{r--} 



{s~ ,5^—0} 



I, (6) 



such that all the amplitudes a^-^z^is'} ^re non-negative ( {r^} 
and {.s^} denote all the possible configurations of t? and s?). 
We can rewrite j^Pf^ ^q) in the following way: 



with 



\Hr'-}) 



2^ Hr^Wl' 



'Nl 



(7) 



E 



=0} 



(8) 



The set {r^}' corresponds to all the configurations of the t[ 



variables such that b 



> 0. Note that each spin state 



IS 



|${t^}) is normalized and X^i^fr-} ^ ^ because 
also normalized. The ground state energy of _ffoff is : 

(*S|i?eff 1*^) = (*^| +BzY ^i^l*') + 

(iJ> i 

(ij) {r'-,r'--}' 

{Hr'llHf^^Mrl). (9) 

is a bounded operator with eigenvalues ±1/2. Therefore 
(${T'^}|i7i^|${T^}) < 1/2 and the equality holds for any 
pair of nearest-neighbor sites (i,j) if |${t^}) = \s^,st = 
N/2) where \s^,st = N /2) is the state with maximum to- 
tal spin (fully polarized) and sf = s^. Hence, expres- 
sion (|9} is minimized for the fully polarized spin configura- 
tion which means that there is family of ground states Y^^) 
that have maximum total spin st and different values of sf,: 



E hr'}\^l 



> (E) 14, ST = N/2). 



(10) 



This proves that there is a fully polarized ferromagnetic 
ground state of -ffoff ■ Note the proof is valid in any dimension 
and is also valid for non-bipartite lattices. When we restrict 
Hcff to the subspace with maximum total spin (st = N/2), 
the operator (Sj • Sj + i) is replaced by 1/2 and the re- 
stricted Hamiltonian becomes exactly the same as the 
one obtained in Ref. 11411 for the strong coupling limit of a 
spinless extended Falicov-Kimball model 12011 . As shown in 
Ref. ClEl, 

the quantum phase diagram of H^^^^ contains a 
ferroelectric phase for < < B'f, i.e., for a non-zero 
value of T-l.. 



The SO(4) symmetry of Hcb implies that the ground state 
degeneracy is higher than the 2st + 1 = + 1 multiplet 
obtained from the global SU(2) spin rotations. Ground states 
with different total spin can also be obtained by making dif- 
ferent global spin rotations for the bands a and b. The spins of 
each individual band will remain fully polarized under these 
transformations, but the relative orientation between spins of 
different bands will change. In particular, the minimum to- 
tal spin will occur when spins in different bands are "anti- 
aligned", i.e., st ~ t^. This implies that the total spin of the 
ground state can take the values: < st < N/2. 



Band a 



Band A 




IPl 

\fPo\ 




mm max 



m„,„, m 



FIG. 1 : Evolution of the magnetization (arrows) on each band under 
the SO(4) transformation U (0). The P vs m plot shows the change 
of the electric dipole moment when the total magnetization evolves 
from the minimum value, obtained for (f) = 7r/2, to the maximum 
value m((j!> = 0) = 1/2. 

To derive the magneto-electric effect, we define \x) as the 
fully polarized ferromagnetic and ferroelectric state obtained 
for B'^^ < Bz < B'f'. Then, the average magnetization and 
electric dipole moment per site, m and p, are given by: 



N 



1 + 



N' 



(11) 



where po = (x|t^|x)/A^. A given ground state in the SO(4) 
multiplet can expressed as = U\x), where U is an el- 
ement of the SO(4) group. In particular, choosing the set of 
transformations C/(0) = e*"^Sj('^j"^'^jti) we get for p(0) = 

M(3u(0)|T-f|5«(0)>/^andTO((/)) = ^ |sTl5n(0))/^: 



m{(f>) 



1 — 4nanh sin + 



N' 



P(0) 



fipo COS ( 



(12) 



where we have used that st = SaT+^bT and that U {(f>) rotates 
the vectors SaT and s^t around the x-axis by angles (f> and —0 
respectively. For the second relation, we have used that: 



U\<t>)T^U (0) =cos0' 



(13) 



and (xK^^lx) = 0- For = 7r/2, the magnetization of the a 
and b bands have opposite sign (see Fig.l) and the total mag- 
netization per site is minimized: nimin = m(7r/2) = \na — 
rib 1/2 (we have taken the thermodynamic limit N oo). The 
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electric polarization can be expressed as a function of m by 
combining Eqs. (I12> : 



v/1 - 4m2 



(14) 



The electric dipole moment is zero for m — rrimin and it in- 
implying that the derivative d\P\/dm 
This important 



creases as y/m — m 
diverges at m = m„iin 



as 1/y'm — m. 

result shows that the interplay between the spin and the or 
bital degrees of freedom can produce an enormous magneto- 
electric effect (see Fig. 1). 
1.0 




2 

z z 

FIG. 2: Zero temperature phase diagram of the two dimensional ver- 
sion of Hcff plus a Zeeman term H^st computed in the spin- wave 
approximation (top) and by exact diagonalization of a 4 x 4 clus- 
ter (bottom). The parameter values are Jo = 0.2J^, — Jy = 
-1.6J°, Ji = -O.IJ,^ and J° = J', + J,/2. 

It is important to verify the stability of the coexisting fer- 
roelectric and ferromagnetic phases away from the limit we 
have considered above. For this purpose, we computed the 
T = diagram of the two dimensional version of i^cfr plus 
a Zeeman term, H^s^, in a spin-wave approximation (top of 
Fig. 2) and by exact diagonalization of a 4 x 4 cluster (bottom 
of Fig. 2). In this case, the values of of Uaa and Ubb are finite 
(Jo = 0.2 J°), and the coexistence of FM and FE obtained for 
= and < < Bf is still present (see Fig. 2). 
For Bz < Bl^ and Hz = 0, the spin-wave phase diagram 
exhibits coexistence of FM and SOO fl4ll . The FM ordering 
is replaced by a spin liquid (short ranged spin-spin correla- 
tions) in the phase diagram computed by exact diagonaliza- 
tion indicating that quantum fluctuations play a major role in 
that regime of parameters. As expected, for large enough Bz 
the system becomes an antiferromagnetic Mott insulator (one 
band is a half-filled and the other one is empty). The exact 
diagonalization shows again that quantum fluctuations gener- 
ate an intermediate phase between the FE-FM state and the 



Mott insulator in which FE and AFM coexist (see bottom of 
Fig. 2). For high values of H^, the system is fully polarized 
and the both T = phase diagrams coincide with the one 
obtained in Ref. |14]. 

The ferromagnetic state can be further stabilized by the in- 
clusion of the ferromagnetic on-site inter-orbital exchange in- 
teraction. For example, the intra-atomic 4/ — 5d exchange 
interaction is about 0.2 eV in EuBg iEbI E4I1 . The other im- 
portant aspect to consider is the role of a finite inter-band hy- 
bridization tab- Exact diagonalization results l23l show that 
the lowest total spin ground state, st ~ ot m = m,nin, is 
the one stabilized after the inclusion of a small tab term. Ac- 
cording to Eq.J14>. this unsaturated ferromagntic state gives 
rise to a divergent magnetoelectric effect. In this situation, a 
small increase in the magnetization produced by an applied 
magnetic field will generate a large increase of the electric 
dipole moment in the way depicted in Fig. 1 . 

In summary, we have shown that the electron-electron 
Coulomb interaction can produce coexisting FM and FE. Both 
phases arise simultaneously from the condensation of excitons 
or particle-hole pairs that exist in two bands with opposite par- 
ity under spatial inversion. The coexistence requires the pres- 
ence of large intra-orbital Coulomb interactions to reduce the 
strength of the antiferromagnetic interaction. We have also 
shown that the coexistence of FE and FM leads to a divergent 
magnetoelectric effect. In the proximity of the ferroelectric- 
ferromagnetic instability, a small magnetic field can produce 
an enormous change in the electric polarization. 
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